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o ■ 

CNJ \ 1. Introduction 

e : 

| ^ ■ In this paper we are interested in the following type of oscillatory integral. Sup- 

. pose S(x, y) is a smooth real-valued function defined in a neighborhood of the origin in 
R 2 , and 4>(x,y) G C^°(R 2 ) is real-valued and supported in a small neighborhood of the 
origin. We define 

U : JsA*)= f e lXS ^4>(x,y)dxdy (1.1) 

_J ' </R 2 

Here A is a real parameter and we want to understand the behavior of Js^(A) as A — > +oo. 
Oscillatory integrals of the form (1.1) and their higher-dimensional analogues come up fre- 
quently in several areas of analysis, including PDE's, mathematical physics, and harmonic 
{Sj . analysis. For example, such oscillatory integrals arise when analyzing the decay of Fourier 
transforms of surface-supported measures such as in [IKM] . We refer to chapter 8 of [S] for 
an overview of such issues. The stability of oscillatory integrals (1.1) under perturbations 
of the phase function S(x, y) is connected to various issues in complex geometry and has 
been considered for example in [PSSt] and [V]. 

O ' 

Since one can always factor out an e lAS, (°> ) ? it does no harm to assume that 
S(0, 0) = 0. If V5'(0, 0) 7^ 0, in a small enough neighborhood of the origin one can 
integrate by parts arbitrarily many times in (1.1) and get that Js^(A) decays faster than 
^ ■ (7/vA _Ar for any N. Hence in this paper we always assume that the origin is a critical point 
H ! for S: that is, we assume that 



£(0,0) = 0, VS(0,0)=0 

In this case, if S(x,y) is real-analytic, using resolution of singularities (see [Gl] for an 
elementary proof) one always has an asymptotic expansion 

oo 

JsM - 5>i^) A_S " + d 'M) ln(A)A" s (1.2) 
j=o 

Here {sj} is an increasing arithmetic progressions of positive rational numbers independent 
of 4> deriving from the resolution of singularities of S. We always assume so is chosen to 
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be minimal such that in any sufficiently small neighborhood U of the origin do(0) or 
d' (4>) is nonzero for some <fi supported in U. In this paper, we will give explicit formulas 
for the leading term of (1.2) once one is in certain coordinate systems which we call 
"superadapted" , in analogy with the adapted coordinate systems of [V]. In the smooth 
case, we will find appropriate weaker analogues. 

Definition 1.1. The oscillatory indexoi S is defined to be sq. If in any small neighborhood 
of the origin there is some (p for which d' ((f)) is nonzero, then we say so has multiplicity 1. 
Otherwise, we say it has multiplicity zero. 

In the case where S is a smooth function whose Hessian determinant at the 
origin is nonvanishing, one can do a smooth coordinate change such that S(x,y) becomes 
±x 2 ± y 2 in the new coordinates. Then one can use the well-known one dimensional 
theory (see Chapter 8 of [S]) and explicitly obtain an asymptotic expansion for Js,<j>(ty', 
the leading term will be given by -^-0(0, 0)A _1 where D denotes the Hessian determinant 

of S at the origin. Hence our concern in this paper will be when D = 0; that is, when S 
has a degenerate critical point at the origin. 

In the real-analytic situation, there is a close relationship between Js,^>(A) and 
the function Is,<t>( € ) defined by 

IsA 6 ) = 4>(x,y)dxdy (1.3) 

J{(x,y):0<S(x,y)<e} 

Analogous to (1.2), when S"(0,0) = the functions Is,<t>( e ) and i_s^(e) have asymptotic 
expansions which we may write as 

oo 

IsAe) ~ j>i(0e r ' + c'M In(e)e^) (1.4a) 

3=0 

oo 

I-sM ~ £(Ci(*)e r ' + C'M) HeV) (1-46) 

j=0 

Analogous to before, {rj} is an increasing arithmetic progression of positive ra- 
tional numbers independent of <p deriving from the resolution of singularities of S, and tq 
is chosen to be minimal such that in any sufficiently small neighborhood U of the origin 
there is some (j> supported in U for which at least one of Co (</>), c' ((f>), Co (</>), or C o (0) is 
nonzero. 

Using well-known methods (see Ch 7 of [AGV]), when the oscillatory index is less 
than 1, corresponding to the degenerate case, r$ = sq and the coefficient of the leading 
term of (1.2) can always be expressed in terms of those of the corresponding terms of (1.4a) 
and (1.46). Often Is,<j> and I-s,<j> are easier to deal with than Js^ due to the absence of 
cancellations which can make it difficult to find lower bounds for Js,^ directly. 
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The Is,4> are closely related to the H$, u defined by 

H s ,u(e) = \{xeU:0< S(x,y) < e}\ (1.5) 

Here U is a small open set containing the origin, and the goal is to understand how 
Hs,u(t) behaves as e — > 0. Many of our results concerning the Is,<f>(t) will immediately 
imply corresponding results about the Hs,u{ e )- Specifically, one chooses <p\ supported in 
U and equal to 1 outside some ^-neighborhood of the boundary of U, then chooses 4> 2 
equal to 1 on U and supported on a 5- neighborhood of U. One compares the theorems 
for Is^iie) and Is,<j> 2 ( e ) an d then lets S go to zero if necessary. In the case where S has 
an isolated zero at the origin, for e small enough the set {x G U : < S(x, y) < e} will 
be a subset of a set where the theorems hold, so one can simply take = 1 and then 
Hs,u(e) = Is,<f( e ) for such e - 

In [V], Varchenko developed some ideas that went a long way towards under- 
standing the case of degenerate real-analytic phase. To describe his work, we need some 
pertinent definitions. 

Definition 1.2. Let S(x, y) = J2 a b s abX a y b denote the Taylor expansion of S(x, y) at the 
origin. Assume there is at least one (a, b) for which s a b is nonzero. For any (a, b) for which 
s a b 7^ 0, let Q a b be the quadrant {(x,y) G R 2 : x > a,y > b}. Then the Newton polygon 
N(S) of S(x, y) is defined to be the convex hull of the union of all Q a b- 

In general, a Newton polygon consists of finitely many (possibly zero) bounded 
edges of negative slope as well as an unbounded vertical ray and an unbounded horizontal 
ray. 

Definition 1.3. The Newton distance d(S) of S(x, y) is defined to be inf{t : (t, t) G N(S)}. 

Throughout this paper, we will use the (ti,t 2 ) coordinates to write equations of 
lines relating to Newton polygons, so as to distinguish from the x-y variables of the domain 
of S(x,y). The line in the plane with equation t\ = t 2 comes up so frequently it has 
its own name: 

Definition 1.4: The bisectrix is the line in the t\-t 2 plane with equation t\ = t 2 . 

A key role in the above theorems as well as our theorems to follow is played by 
the following polynomials. 

Definition 1.5. Suppose e is a compact edge of N(S). Define S e (x,y) by S e (x,y) = 
S( a b)ee s abX a y b - In other words S e (x,y) is the sum of the terms of the Taylor expansion 
of S corresponding to (a, b) G e. If S(x,y) is real-analytic, we use the same terminology 
when e is the vertical or horizontal ray of N(S). 

In [V], Varchenko showed that when S is real-analytic the oscillatory index is 
always at most -rrgr, and that there is necessarily a coordinate system in which it is 
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actually equal to He also showed that the coordinate change to such coordinates can 
always be made of the form (x, y) — ■> (x, y — f(x)) or (x, y) — > (x — f(y),y) for real analytic 
/. Coordinate systems where d(S) achieves the maximum possible value are referred to 
as "adapted coordinates" . He also showed that the multiplicity of the oscillatory index is 
equal to 1 if and only if there are adapted coordinates where the bisectrix intersects N(S) 
at a vertex. Otherwise the multiplicity is 0; the leading term of (1.2) will not have the 
ln(A) factor in it. The issue of finding an expression for the leading coefficient do ((f)) or 
d' (<p) is not treated in [V]. However, in the case where the S(x,y) has a critical point of 
finite Milnor number at the origin, it is shown in [V] that the leading coefficient do((f>) or 
d'o((f>) is some fixed multiple of 0(0, 0) depending on the phase; precisely which multiple is 
not determined. These results were later extended in [Sh] to convex finite-type functions. 

Smooth analogues are proven in [IM] and [IKM] . In [IM] it is shown that adapted 

coordinates exist in the smooth case. In [IKM] it is shown that in smooth adapted coordi- 

_ i 

nates one has the estimate | Js,<t>(*)\ < Cln |A||A| d( - s ^ for large |A|. There are also operator 
versions of such results. For example, in [PS] it is proven that 1 1 J e lXS ( x,y > (f)(x, y)f(y) dy\ \ L * 

< C\\\~™&\\f\\ L 2 for real analytic phase, where S(x,y) = S(x, y) - 5(0, y) - S(x, 0). The 
exponent is sharp. Generalizations to smooth phase were proven in [R] and [G2]. 

These results use subdivisions into curved regions as will be done here. However, there are 
significant differences since one gets stronger results for operators; in particular, adapted 
coordinates are not needed. 

Our theorems below will require us to be in certain adapted coordinate systems 
which we call "superadapted" coordinate systems: 

Definition 1.6. One is in superadapted coordinates if whenever e is a compact edge of 
N(S) intersecting the bisectrix, both of the functions S e (l,y) and S e (—l,y) have no real 
zero of order d(S) or greater other than possibly y = 0. 

It can be shown that an equivalent definition is obtained by stipulating the same 
condition on S e (x,l) and S e (x,—1) instead of S e (l,y) and S e (—l,y); we choose the y- 
variable for definiteness. In section 7 we will prove any phase function can be put in 
superadapted coordinates using some ideas from two-dimensional resolution of singulari- 
ties. 

Lemma 1.0. In a superadapted coordinate system, a critical point of S(x, y) at the origin 
is nondegenerate if and only if d(S) = 1. 

Proof. Write S(x,y) = ax 2 + bxy + cy 2 + 0(|x| 3 + \y\ 3 ). The only way d(S) could be 
greater than 1 is for either a and b to both be zero, or for c and b to both be zero. In 
either case, the Hessian at the origin is zero and the phase is degenerate. So we assume 
that d(S) = 1, and we will show that in a superadapted coordinate system the phase is 
nondegenerate. 
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First consider the case where N(S) has a vertex at (1, 1). Then 6 ^ and either 
a or c is zero. Suppose a = but c 7^ 0; we claim that this implies the coordinate system 
is not superadapted. For in this case there is an edge e connecting (1, 1) and (0, 2). Then 
S e {x,y) = bxy + cy 2 and thus S e (l,y) = by + cy 2 has a zero at — \ 7^ 0, inconsistent 
with the definition of superadapted. Thus in a superadapted coordinate system, if N(S) 
has a vertex at (1,1) and a = 0, then c = and thus the Hessian is nonzero at the 
origin; the phase is nondegenerate. The case where c = but a ^ leads to a similar 
contradiction. We conclude that if N(S) has a vertex at (1, 1) and d(S) = 1 then the phase 
is nondegenerate. 

Next, consider the case where (1, 1) is in the interior of an edge e of N(S). In this 
case, the endpoints of e are (2, 0) and (0, 2). Hence a and c are nonzero. In a superadapted 
coordinate system, one must have that S(l,y) = a + by + cy 2 has no real zeroes other than 
y = 0. Since o ^ 0, this is equivalent to a + by + cy 2 having no real zeroes at all, which 
happens exactly when b 2 < Aac. This is equivalent to the Hessian determinant at the origin 
being nonzero, and thus the phase is nondegenerate in this situation too. This completes 
the proof of Lemma 1.0. 

We now come to our theorems. We will use the shorthand d to denote the 
Newton distance d(S). If N(S) intersects the bisectrix in the interior of an edge, bounded 
or unbounded, we denote this edge by eo and its slope by — — , where < m < 00. We use 
the shorthand So(x, y) to denote S eo (x,y). 

Theorem 1.1 is our main result. It gives explicit expressions for the leading term of 
(1.4a) in superadapted coordinates. Applying the theorem to — S gives analogous formulas 
for the expansion (1.46). As indicated above, (1.4a) — (1.46) directly imply formulas for 
the leading term of the asymptotic expansion (1.2) in the degenerate case; these are given 
in Theorem 1.2. 

Theorem 1.1. Suppose S(x,y) is a smooth phase function in superadapted coordinates 
with d > 1. If the function 0(x, y) is supported in a sufficiently small neighborhood of the 
origin, then r$ = ^ and the following hold. 

a) Suppose the bisectrix intersects N(S) in the interior of a compact edge. Define the 
function SQ(x,y)~d to be S (x,y)~d when S (x,y) > and zero otherwise. Then we 
have 

lim^4^ = (m+l)-V(0,0) f°° (St(l,y)-* + S+(-l,y)-*)dy (1.6) 

In particular, if S(x,y) is real-analytic then the coefficient c' ((f)) in (1.4a) is always zero 
and cq{4>) is given by (1.6). 

b) Suppose the bisectrix intersects N(S) at a vertex (d, d). Let SddX d y d denote the corre- 
sponding term of the Taylor expansion of S; hence sad 7^ 0. Denote the slopes of the two 
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edges of N(S) meeting at (d, d) by s\ and s 2 , where — oo < s 2 < s\ < 0. Then 

lim ^{ = V (S)\s dd \-h(0M-^-r ~ -^r) (1-7) 
e~d ln(e) si — 1 S2 — 1 

Here 77(5") = 4 if s dd > and d is even, rj(S) = 2 if d is odd, and 77(6') = if s dd < and d 
is even. 

c) Suppose S(x,y) is real-analytic and the bisectrix intersects N(S) in the interior of the 
horizontal ray. Write So(x,y) = a(x)y d where a(x) is real-analytic. Let a(x) denote the 
one-dimensional measure of {y : < a(x)y d < 1}. In particular, a(x) = \a(x)\~d when d 
is odd. Then c' (<j)) = 0, and cq{4>) is given by 

/oo 
a(x)(j)(x,0)dx (1.8) 
-00 

The case where the bisectrix intersects the interior of the vertical ray has the analogous 
formula. 

As can be seen, parts a) , b) and c) give quite different formulas. Correspondingly, 
in our subsequent theorems we break up into three cases. Case 1 is when the bisectrix 
intersects N(S) in the interior of a compact edge, case 2 is when the bisectrix intersects 
N(S) at a vertex (d, d), and case 3 is when the bisectrix intersects N(S) in the interior of 
one of the unbounded rays. Notice that in cases 1 and 3, for a given (f>(x, y) the expressions 
of Theorem 1.1 depend only on Sq(x, y), and that in case 2 it depends on SddX d y d as well 
as the slopes of the edges of N(S) meeting at (d, d). 

Our next theorem gives the oscillatory integral version of Theorem 1.1 for the 
real-analytic case. 

Theorem 1.2. Assume S(x, y) is real-analytic and is in superadapted coordinates with 
d > 1. Then so = \. In case 1 and 3, the coefficient d' (4>) of (1.2) is always zero and do(</>) 
is given by 

d o (0) = ^(e^c o (0) + e-^C o (0)) (1.9a) 



In case 2, one has 



d'M = -^-(e^c'M + e-^C'M) (1-96) 



Proof. We will be sketchy here since the method for proving Theorem 1.2 from Theorem 
1.1 is well-known; we refer to chapter 7 of [AGV] for more details. Note that 

/•OO poo 

JsAV= (dJsAt)y X£ de+ {dj. s ^e))e-^de (1.10a) 
Jo Jo 
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By [F], if 7 G C C (R) with j(t) = 1 near and if a > — 1, for any / we have 

e iX ^Htri(t)dt=^ I ^L + 0(X- 1 ) (1.106) 

Inserting (1.4a) and (1.46) into (1.10a) and using Theorem 1.1 and (1.10b) gives the the- 
orem. 

Comment 1. One does need that d > 1 for Theorem 1.2 to hold. When S(x,y) has a 
nondegenerate saddle critical point, there are coordinates where S(x, y) = xy. This falls 
under case 2, and one has that Cq(0) = Cq(0). This means that the two terms of (1.96) will 
cancel. And in fact when 0(0, 0) 7^0, \IsA e )\' \I-sA € )\ ~ l eni ( e )|j while | Js,</>(A)| ~ A -1 . 

Comment 2. In cases 1 and 2, the expressions (1.9a) and (1.96) for d o (0) and d' ((j)) will 
always be nonzero when d > 1 and 0(0, 0) ^ 0. This is because the expressions given by 
Theorem 1.1 for co(0), Co(0), Cq(0), and Cq(0) are real multiples of 0(0, 0), while the ratio 
of the e l 2d and e - *^ factors is never real when d > 1. 

Comment 3. In any dimension, when the phase satisfies an appropriate nondegeneracy 
condition there are reasonably explicit formulas for the leading coefficient of the leading 
term of the asymptotic expansion of oscillatory integrals such as (1.1). Such formulas are 
proven in [DS] and [DNS]. 

Next, we give some less precise C°° analogues for the IsA e )- The lower bounds involve 

I\S\A*) = ^At) + I-sA e ) = I{( x ,y):\S(x,y)\<e} ^ V) dx d V ' 

Theorem 1.3a. Suppose now that S(x,y) is a smooth phase function in superadapted 
coordinates with d > 1. If is supported in a sufficiently small neighborhood of the origin, 
then there is a positive Bs,4> such that: 

In cases 1 and 3 one has 

IsAz) < B s ,*e* (1.11) 

In case 2 one has 

IsAe) < B sA Hz)\t~ d (1-12) 

One has some analogous lower bounds for I\s\A e )- They are sharp in cases 1 and 2, and 
almost sharp in case 3. (Sharp lower bounds do not hold in general in case 3, as explicit 
examples show). 

Theorem 1.3b. Suppose we are in the setting of Theorem 1.3a). Suppose also that 
0(0,0)^0. 

In case 1 there exists a Ag^ > such that for e sufficiently small we have 

I\s\A e ) > A s,^ (1.13a) 
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In case 2, one similarly has 



I\s\A € ) > ^l ln ( e )l e " 



(1.136) 



In case 3, one has analogous almost-sharp lower bounds, at least if 4>(x,y) is nonnegative. 
Namely, for any 5 > one has 



The next lemma will be used to show that the three cases of superadapted coor- 
dinates are mutually exclusive. 

Lemma 1.4. Assume S(x, y) is smooth and is in case 3 of superadapted coordinates with 
d > 1. Then for any M one can find a smooth function Sm(x, y) such that Sm — S has a 
zero of order at least M at the origin, but such that in a small enough neighborhood U of 
the origin one has 



Proof. Suppose for example that the bisectrix intersects N(S) in the interior of the 
horizontal ray. Then Sm(x, y) = S(x, y) +x M agrees with S(x, y) to order M at the origin 
and has Newton distance less than that of S(x,y). The Newton distance is also greater 
than 1 for M large and the relevant polynomials (S , M)e(l, y) and (5 , M)e(— 1, y) have zeroes 
of order at most 1. Hence SM{x,y) is in case 1 superadapted coordinates and one can 
apply Theorem 1.1a to conclude that (1.14) is finite for a small enough neighborhood U 
of the origin. 

Lemma 1.5. Any smooth degenerate phase S(x, y) can be put in superadapted coordi- 
nates in exactly one of cases 1, 2, or 3. 

Proof: In section 7 we will show that one can always put S(x,y) into some superadapted 
coodinate system. Equations (1.11) (for S(x,y) and —S(x,y)) and (1.136) cannot simul- 
taneously hold, so a case 2 coordinate system cannot be put in a case 1 or 3 coordinate 
system. Suppose S(x, y) has a case 3 coordinate system as well as a case 1 coordinate 
system; we will derive a contradiction. Since it has a case 3 coordinate system, we may 
adjust S(x,y) to arbitrarily high order and cause (1.14) to hold. In its case 1 coordinates, 
an adjustment of high enough order will not affect the fact that it is in case 1 and thus 
equation (1.11) will still hold. Using the well-known relationship between LP norms and 
distribution functions (on the function [I\s\ ,<^(e)) _1 ), one then gets that the integral (1.14) 
in the new coordinates is infinite, a contradiction. Thus the three cases are mutually 
exclusive. 

For oscillatory integrals with smooth phase, one has some analogues of Theorem 1.3. It 
should be noted that Theorem 1.6a can be proved using the results of [IKM]. 



I \s\A e ) > A S,4>,8£ 



(1.13c) 




(1.14) 
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Theorem 1.6a. Suppose S(x, y) is smooth and is in superadapted coordinates with d > 1 
In cases 1 and 3 as A — * oo one has 



|Js,*(A)| <C\—* 



In case 2 one has 



\JsA*)\ <c\—*h\) 

Theorem 1.6b. Suppose y) is nonnegative with 0(0, 0) > 0. 
In case 1, one has 



In case 2, one has 



In case 3, for any 5 > one has 



lim sup 

A^oo 



lim sup 

A^oo 



A d 



A-aln(A) 



> 



> 



lim sup 

A^oo 



A d" 



= oo 



(1.15a) 
(1.156) 



(1.16a) 



(1.166) 



(1.16c) 



Although we will not prove it here, it can be shown with some additional argument that 
the conditions on <j>(x,y) in (1.16a) — (1.166) can be weakened to just that 0(0,0) ^ 0. 

In all of cases 1, 2, and 3 we will divide the domain of integration of the expressions 
(1.1) and (1.3) for Jg^ and Is,</> into 4 parts, depending on whether or not x and y are 
positive or negative. Adding the resulting formulas and estimates will give the theorems. 
Without loss of generality we will always focus on the x, y > as the other quadrants are 
always dealt with the same way. Hence our goal is to understand Ig ^ and Jg^, where 

J s^( e ) = / <f){x,y)dxdy (1.17a) 

J{(x,y):x>0,y>0, 0<S(x,y)<e} 

44 X )= I e lXS ^<P{x,y)dxdy (1.176) 

J{(x,y):x>0,y>0} 

In turn, the domains of (1.17a) — (1.176) will be written as the union of various "curved 
triangles" (such as those of Lemma 2.0 below). On a given curved triangle, one typically 
Taylor expands S(x, y) or one of its derivatives and then uses Van der Corput-type lemmas 
in the x or y direction to get a desired estimate. For the oscillatory integrals, the traditional 
van der Corput (see Ch 8 of [S]) is used, while for sublevel integrals the version of [C] is used. 
Van der Corput-type lemmas have been considered in some detail, such as in [ArKaCu] 
and [CaCW] , as well as the early work of Vinogradov [Vi] . We refer to [CaCW] for further 
results and references. 
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Throughout this paper, we will often have a constant C appearing on the right- 
hand side of an inequality. This always denotes a constant depending on S and 0. Occa- 
sionally we will need further constants C, C", etc which also depend on S and 4>. 

2. Some useful lemmas for Cases 1 and 2. 

Suppose G is an open subset of R 2 . Then throughout the course of this paper 
we will make frequent use of Ig ^(e) and Jg ^(A) defined by 

IsAe) = I <P(x,y)dxdy (2.1a) 

J{(x,y)eG:0<S(x,y)<e} 

Js,4 X )= f e* xs ^<P(x,y)dxdy (2.1b) 
Jg 

A certain type of G comes up in several contexts in this paper, and relevant estimates we 
need for Ig ^(e) and Jg^W are given by the following lemma. 

Lemma 2.0. Suppose for some A, m > 0, < 8 < 1, we let G = {(x,y) G [0,5] x [0,6] : 
< y < Ax m }, and suppose there are nonnegative integers a and b with a > b and a > 2 
such that for some constant Cq the following holds on G. 

d b y S(x,y) >C x a (2.2a) 

If b = 1, assume also that 

d 2 y S(x,y) <C x a - m (2.26) 

If b = 0, instead of (2.26) assume also that a > m + 1 and that for some constant C\ we 
have 

d x S(x,y) > Cix a ~\ d 2 x S(x,y) < C x a ~ 2 (2.2c) 

Then for some ( a b > and C depending on S, 0, and Cq (and C\ if 6 = 0), if the support 
of cj) is contained in [—6, 5] x [—5, 5] one has 

< Ce^A Cab (2.3a) 

\J^(\)\<C\\\-^A^» (2.36) 

Proof. We first consider the case where 6 > 0; we will do the 6 = argument afterwards. 
By (2.2a) and the Van der Corput lemma in the y direction (see [C] for example) one has 
that for a given x we have 

|{y:|S(s,y)|<e}|<Ce*a;-t (2.4) 
As a result, if G x denotes the vertical cross section of G at x, of length Ax m , then one has 

\{y E G x : \S(x,y)\ < e}\ < Cmm(Ax m , e^x~^) (2.5) 
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Consequently, we have 

rS 

\ I s<f>\< C mm(Ax m ,e^x~^)dx (2.6) 
Jo 

It is natural to break the integral (2.6) into two parts, depending on whether or not 

Ax m > e^x~%. The two quantities are equal at xo = A~ a + mb e a + mb . The left integral 
becomes 

f ° A 1 a-b m+1 

/ Ax m dx = = i«+»- k e»+'»' (2.7) 

Jo m + 1 m + 1 v ' 

The right integral is computed to be 



i i f°° 

e b / x~ dx < e b / x - ' 5 ' da; 



6 1 6 , a-b m+1 

-e&Xg = -i»+-»»e»+" 1 ' (2.8) 



a — o a 



Adding together, we obtain that |/g^| < CM. a + mb e a + m! > as needed. 

The estimates for «/g^(A) for 6 > 2 are done in a similar fashion. First suppose 
6 > 2. Then one can use the usual Van der Corput lemma (see [S] ch 8) in the y direction 
to obtain 

e iXS M(j)(x,y)dy\<C\X\-^x-^ (2.9) 



L 



G x 

-1 



This is the analogue to (2.4) with e replaced by |A| . As a result, similar to (2.6) we get 

\J^\<C mm(Ax m ,\X\-^x-^)dx (2.10) 
Jo 

The result is 

f i a — b m-\-l 

\Js^\ < CA^+^\X\~^+^ (2.11) 

This gives (2.36). We next prove (2.36) when 6=1. If one integrates by parts in the y 
variable one gets several terms each of which can be bounded using (2.2a) and (2.26). If 
one works it out, one gets that these terms are bounded by C|A|- 1 x-^ a+m) . It is thus 
natural to split the integral into two parts at the xq satisfying |A| _1 xo a+m = 1, in other 
words, at Xq = \ X\~ a + m . We accordingly write G = G\ U (?2, with G\ the portion where 
xo < |A|~ a + m . We then have 



m + 1 



Js^X) KCldl <C'A\X\-— (2.12) 

As for jf^(A), we integrate by parts in y, obtaining the C\X\~ 1 x~^ a+m ' ) factor, and one 
gets that 

\Js,%W\<cf lAp 1 *-^) 
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1 r Ax 



<C I iX^x-^+^dydx 

J\\\~a JO 

= CA I , Wl^x^dx 



|,\| a+m 

<CA\X\-^ (2.13) 
Adding (2.12) to (2.13) give the oscillatory integral estimates for 6=1. 

We now consider the oscillatory integral when 6 = 0. Here we do the integrations 
by parts in the x direction. This time, by (2.2a) and (2.2c) an integration by parts incurs 
a factor of |A| _1 x _a . Hence we subdivide G = Gi U G 2 , where G\ = {(x,y) G G : < x < 
|A| _ «}. Note that the measure of G\ is ^lAI -12 ^ 1 , so that 



J^{\)\<CA\\\-^ (2.14) 



m + 1 



For the G 2 piece one obtains 



\4U X )\ <C [ \M~ lx ~ a = f 1 I"" IM'^dydx 
Jg J\x\-a Jo 

1 f , 

J\X\-a 



= cwr 1 / . Ax m - a dx 



m+1 



= CA\X\-— (2.15) 

(For the last equality we use the hypothesis that a > m+1). Adding (2.14) to (2.15) gives 
the estimate we seek, (2.36). Lastly, we prove the bounds for Ig ^(e) when 6 = 0. In this 
case, since \S(x,y)\ > Cox a , we have 

1^1 < C\{(x,y) E G : C x a < e}\ < C'{(x,y) E G : x < e^} = C" Ae^ (2.16) 

This concludes the proof of Lemma 2.0. 

In Lemmas 2.1 and 2.2 below, S(x, y) is a smooth phase function in Case 1 or 2 of 
superadapted coordinates with d > 1. If there is a compact edge E of N(S) such that the 
bisectrix contains either the upper vertex of E or an interior point of E, then we denote 
the equation of this edge by t\ + mt2 = a, and for some large but fixed number we let 
A\ = {(x, y) £ [0, 1] x [0, 1] : y < jjX m }. Similarly, if there is some edge E' with equation 
t\ + m't2 = a' such that the bisectrix contains either the lower vertex of E' or an interior 
point of E', we let A 2 = {(x,y) G [0,1] x [0, 1] : x < -^y~}- Note that in case 1 both 
A\ and A 2 exist. We focus our attention on I^i 6 ) an d <^s^(A). The relevant information 
about them (if they exist) is provided by the following lemma. 
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Lemma 2.1. There exists an 77 > such that if the support of y) is sufficiently small, 
then for z = 1 , 2 we have 

|/^(£)|<Ce3JV-" (2.17) 

\J&{\)\<C\\\-tN-'> (2.18) 

Proof. By symmetry, it suffices to prove the bounds for A\. Let (a, b) denote the lowest 
vertex of E. Thus b < a. Since (a, 6) and (d, d) are both on E, we have a + mb = (1 + m)d 
or i = ■ We will show that the hypotheses of Lemma 2.0 hold for these values of 
a, 6, and m, setting A = A -1 . Since \ = rn . + \ , Lemma 2.1 will follow. For a large but 

' ' jo _ d a+mb ' ° 

fixed M, we write S(x,y) = J2 p <m q <M s pq xP V q + Em(x,v)- By standard estimates, for 
< a, (3 < M we have 

\d^E M (x,y)\ < C(\x\ M - a + \y\ M ^) (2.19) 

We can write 

d b y S(x,y)= s' pq x p y q + d b E M (x,y) (2.20) 

p<M, q<M-b 

Here s' a0 7^ 0. We next show that the sum in (2.20) is dominated by the term s' a0 x a in a 
sufficiently small neighborhood of the origin. To this end, note that (a, 0) is a vertex of 
the Newton polygon of dyS(x,y), and that a horizontal ray and an edge of this Newton 
polygon with equation t± + mti = a intersect at (a, 0). As a result, for any term s' pq x p y q 
in the sum of (2.20) other than s' a0 x a , either p > a, or p < a, q > 0, and p + mq > a. 
Correspondingly, we let Ti = {(p, q) : p < a, < q < M — b, p + mq > a} and T2 = 
{(p, q) : a < p < M, < q < b — a, (p, q) 7^ (a, 0)} and we rewrite (2.20) as 

d b y S(x, y) = s' a0 x a + J2 VV + E V P y q + d v E M(x, V) (2.21) 

We examine a given term s' pq x p y q in the T\ sum. Since (x, y) is in the domain A±, one 
has y < As a result, Is^a^y 9 ! < |Sp (J A^~' 7 x p+m ' 3, | < \s' pq N~ q x a \. Thus so long as is 
chosen sufficiently large (depending on M) , we may assume that the absolute value of the 
T\ sum is at most ||Sa l xa - We next examine a term s' pq x p y q in the T 2 sum. Here we have 
\ s ' pq xP y q \ < kpg^KM + \y\)- Hence if we are in a sufficiently small neighborhood of the 
origin, we can assume the absolute value of the T 2 sum is also at most ||s' a ol^ a - Lastly, 
by (2.19), in a sufficiently small neighborhood of the origin \d y EM{x,y)\ is also bounded 
by tjlsaol^"- Putting these together, we conclude that on the domain A\ we have 

\d b y S(x,y)\>^\s' a0 \x a (2.22) 

Note that when b > 1, (2.22) ensures hypotheses of Lemma 2.0 hold. Thus one may apply 
Lemma 2.0, giving Lemma 2.1. When 6 = 0, in order to apply Lemma 2.0 one needs also 
that a > m + 1 and that (2.2c) holds. But taking an x derivative of S(x,y) just shifts 
the Newton polygon to the left by 1, so exactly as in (2.22) the first and second derivative 
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conditions of (2.2c) will hold. As for the requirement that a > m + 1, note that (d, d) and 
(a, 0) are on the line ti +mt2 = a and that d > 1. This means a = a+mO = d+md > 1+m 
as needed. 

Lastly, we show that the supplemental hypothesis (2.26) holds when 6=1. It 
helps to view things in (x,y') coordinates where (x,y) = (x,x rn y / ). The condition (2.26) 
becomes that \dy,S(x, x m y')\ < Cx a+m . Also, since the terms s pq x p y q of S(x,yYs Taylor 
expansion with minimal p + mq (= a) are exactly the terms of SE{x,y), the expansion 
S(x,y) = J2 p <m q <M s pq xP U 9 + EM^x.y) becomes of the following form, where Tm is a 
polynomial in y' and a fractional power of x. 

S(x, x m y') = x a S E (l, y') + x a+e T M (x, y') + E M (x, x m y') (2.23) 

Using (2.23) and the error estimates (2.19) we have \dy,S(x, x m y') \ < Cx a . But (a, 6) = 
(a, 1) is on the line ti +mt2 = a and therefore a + m = a. This gives the second derivative 
bounds of (2.26), and thus the hypotheses of Lemma 2.0 hold. This concludes the proof of 
Lemma 2.1. 

Lemma 2.2. Let S(x, y) be a smooth phase function in case 1 or 2 superadapted coor- 
dinates with d > 1. Suppose e is a compact edge of N(S) intersecting the bisectrix and 
has equation given by t\ + mt 2 = a. Define B = {(x,y) e [0, 1] x [0, 1] : ^ < y < Nx m }, 
where is some large but fixed constant. Then if the support of (p(x, y) is sufficiently 
small, depending on A, one has the estimates l-ff^e)! < C^e^ and | J^(A)| < CV|A|~d. 

Proof. In the above (x,y f ) coordinates one has 

/| (e)= / x m <P(x,x m y')dxdy' 

J{(x,y')e[0,l]x[N- 1 ,N]:0<S(x,x m y')<e} 

J|^(A)= I e lXS{x ' xmy ' ) x m <P(x,x m y')dxdy' 

J[0,l]x[N-\N] 

In view of (2.23), the zeroes of S e (l,y') might be expected to play a significant role in the 
analysis. To this end, we assume that A is large enough so that any zeroes of S e (l,y') for 
y' > are in [A -1 , A], and denote these zeroes by z\, Zk (if there are any). Let vi, ...Vk 
denote the orders of these zeroes, and let U denote the interval [zi — ^, Zi + ^]. By (2.23) 
and the error term derivative estimates (2.19) we can assume that on the (sufficiently 
small, depending on A) support of 4>(x,x m y'), if (x,y') G [0, 1] x Ii then 

\dpS(x,x m y')\ > Cx a (2.24a) 

In the case that vi = 1, we similarly have 

\dlS{x,x m y')\ < Cx a (2.246) 
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We now translate this into the orginal (x, y) coordinates. Each [0, 1] x Jj becomes a set 
Di of the form {(x, y) E [0, 1] x [0, 1] : {z % - ±)x m < y < (z l + ^)x m , and on D t (2.24a) 
becomes 

\d£S(x,y)\ > Cx a ~ mVi (2.25a) 
In the case of Vi = 1, (2.246) becomes 

\dyS(x,y)\ < Cx a - 2m (2.256) 

Since we are superadapted coordinates, < Vi < d. Thus (m + l)vi < (m+l)d. Since (d, d) 
is on the edge e, we have a = (m + l)d. Thus (m + l)vi < a or a — mvi > Vi. As a result, 
the sets Di have vertical cross sections of length j^x m and satisfying (2.25a) — (2.256) with 
a — mvi > Vi. Hence after doing a coordinate change of the form (x, y) — > (x,y — f(x)), 
we are in the set-up of Lemma 2.0 and for some i] > we get 

1^1 < Ce^N~ v (2.26a) 



Js,4>\ < C\X\-^N~ V (2.266) 



Since = \ the above becomes 

a d 



\I^\^ Ce ~ dN ~ V ( 2 - 27a ) 
\Js,4,\<C\X\~^N- v (2.276) 

Next, write [iV -1 ,iV] — Uj/j as the union of intervals Jj. Then since S e (l,y') has no 
zeroes on any Jj, by the expansion (2.23) and the error derivative bounds (2.19), if 5n is 
sufficiently small then on [0, 8n] x Jj we have 

\d x S(x,x m y') \ > Cnx*' 1 \d 2 x S{x,x m y')\ < C' N x a ~ 2 (2.28) 

Separating at x = |A|~« and integrating the right portion by parts in x using (2.28) as in 
the proof of Lemma 2.0 gives 



f e lXSi - x ^ my ' ) x m (t){x } x rn y')\dxdy' <C^\X\-^ (2.29a) 

J[0,S N ]xJi 



Converting back into (x, y) coordinates and using that \ = (2.29a) becomes the 

following, where Ei denotes the set [0, 5n] x Jj in the (x, y) coordinates. 

| f e lXS ( x ^<J)(x,y)dxdy\<C^\\\-^ (2.296) 

JEi 

Moving now to /^^(e), by (2.23) one has \S(x, x m y') \ > Cxx a on [0, 8n] x Ji for sufficiently 
small 5n > and some CV (not necessarily the same constant as above). As a result, 

x m (j)(x,x rn y') dx dy'\ 

{{x,y')e[0,S N }xJ x :0<S(x,x™y')<e} 
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< I / x m (f)(x,x m y')dxdy'\ 

JUx,y')e\0,S N ]xJ i :0<C N x°'<e} 

'AT / 
JO 



'{(x,y')e[0,5 N ]xJi:0<C N x<*<£} 
/•e a 

< C' M / x m dx 



C^e^ (2.30a) 



Again going back to the (x, y) coordinates and using that \ = 222 ^~; we conclude that 

|/ <p(x,y)dxdy\ < C N e^ (2.306) 

J{(x,y)eE,:0<S(x,y)<e} 

Lemma 2.2 now follows by adding (2.296) - (2.306) to (2.27a) - (2.276). 



3. Case 1 proofs. 



Assume now that we are in Case 1. We start by proving the upper bounds for 
smooth phase. 

Theorem 3.1. The right-hand sides of (1.11) and (1.15a) hold. 

Proof. In case 1 of superadapted coordinates, the domain of integration of /^(e) or 
Jg^A) is the union of Ai, A2, and B : where A\ and A2 are as in Lemma 2.1 and where B 
is as in Lemma 2.2 for the edge of N(S) intersecting the bisectrix. Thus the theorem follows 
by fixing some N and adding the inequalities of Lemmas 2.1 and 2.2 to the corresponding 
inequalities for the other quadrants. 

The next lemma will be useful in getting the formulas for the real-analytic case. 

Lemma 3.2. Suppose S(x,y) is a smooth case 1 phase function like before. Then there 
is a natural number D < d and a neighborhood U of the origin such that if <f>(x, y) is 
supported in U and is zero on a neighborhood of the origin, then |is,^(e)| < Ce& . 

Proof: It suffices to fix some iV and show that each -fg^(e), -^^(e), -^s^( e ) satisfies 
the upper bounds, where the Di and correspond to the edge of N(S) intersecting the 
bisectrix. We start with Ig^. Without loss of generality we may take i = 1. Since <f)(x,y) 
is zero on a neighborhood of the origin, there is some 5 > such that <fi(x,y) is zero for 
(x, y) e Ax with < x < 5. (2.22) says that dyS(x,y) is bounded below on A±, where 
e < d denotes the ^-coordinate of the lower vertex of the edge of N(S) intersecting the 
bisectrix. If e > 0, by the Van der Corput lemma in the y direction, we have 

\{y e A\ : < S(x, y) < e}| < Ce= (3.1) 

Thus we have 

|/^(e)| < C f \{y e A\ : < S(x,y) < e}\dx < Ct^ (3.2) 
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This is the desired estimate for e > 0. If e = 0, then (2.22) says that S(x, y) is bounded 
below on the support of the integrand of Ig 1 ^) and then (3.2) holds trivially. Thus we 

have the desired bounds for the Ig 1 ^). The Ig^ie) are dealt with in a similar way. This 
time, one uses (2.24a) to obtain |i^(e)| < Ce v i . 

Lastly, we look at the |/J (e)|. As mentioned below (2.296), by (2.23) one has 
\S(x,y)\ > Cx a on each E^. Hence since (f>(x,y) is zero on a neighborhood of the origin, 
\S(x, y)\ is bounded below on the support of the integrand of Is%( e ) and (3.2) again holds 
trivially. This completes the proof of Lemma 3.2. 



We now move to the case 1 formulas of Theorems 1.1. We write = 0i + 02, 
where 0i is supported in a smaller neighborhood of the origin and 02 is zero on a neigh- 
borhood of the origin. By Lemma 2.3, lim e ^o Js, *j^ = and therefore lim e ^o ^'l^ = 

lim e ^o Js '^i ■ Hence when proving the limit (1.6) one can always replace by 0i at will, 

e d 

regardless of how small the support of 0i is. 

Our strategy will involve fixing iV and analyzing the Ig^, where 0i has small 
support depending on N. Lemma 2.1 and (2.27a) will ensure that the contributions of the 
Igfa and Ig^ will be 0(N~ ri ) smaller than that of the Ig^ as e goes to zero. Letting 

N go to infinity will give lim e ^ -^r — = Si nm e^o -^r — • Adding up the latter limits 

e d e d 

along with their analogues in the other three quadrants will give (1.6). 

Proof of (1.6). It suffices to assume 0(0,0) > as the general case can be obtained 
by writing = 0' — 0" where 0'(O,O) and 0"(O,O) are positive. We work in the (x,y') 
coordinates like above. Namely, we write 

I%{e)= f x m cf>(x,x m y')dxdy' (3.3) 

J {(x,y')e[0,l]x Ji:0<S(x,x m y')<e} 

Ji was defined so that So(l, y') has no zeroes on Jj. As a result, by (2.23) (and using (2.19) 
to deal with E M {x,x m y')), we may let 5 > such that on [0,5] x Jj either S(x,x m y') is 
negative, or we have 

< (1 - ±)S(x, x m y') < x«S (l, y') < (1 + ±)S(x, x m y') (3.4) 
Shrinking 5 further if necessary, we assume 5 is small enough so that on [0,6] x J i we have 

<f>(x,x m y')< (1 + 1)0(0,0) (3.5) 



As described above, one can multiply 4>(x,y) by a cutoff function supported on \x\ < 5 

I Ei (e) 

without affecting lim e ^ -^r — • Hence we assume <fi(x,y) is supported on \x\ < 5, and 
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that the multiplying cutoff was chosen so that for (x, y') E [0, 5] x J i? equation (3.5) holds. 
We also assume the multiplying cutoff was chosen so that for (x, y') E [0, |] x Jj one has 

(l-^)<P(0,0)< ( j ) (x,x m y') (3.6) 

We now proceed to our main estimates. If S(x,x m y') is negative on [0,5] x J i? Ig^e) 
becomes zero. If on the other hand (3.4) holds, then by (3.4) and (3.5) one has 

I^(e) < (1 + -^(0,0) / x m dxdy' (3.7a) 

iV J{(x,y>)e[0,S]xJ l :0<x°<So(l,y')<(l + -k)e} 

On the other hand, by (3.4) and (3.6) we also have 

(e) > (1-1)^(0,0) / x m dxdy' (3.76) 

We now change coordinates from x to x' = x m+1 in the integrals (3.7a) — (3.7b). Using 
the fact that = d, which follows from the fact that (d, d) is on the line ti + mt 2 = a, 
(3.7a) — (3.76) become the following, where 5' = 5 m+1 . 

/|Ue)<(m+l)- 1 (l + -|-)0(O,O) [ dxdy' (3.8a) 

I^yim + iy^l-hmO) [ dxdy' (3.86) 

Doing the x integrals first, (3.8a) — (3.86) become 

/|;(e)<(m+ 1)^1 + ^)0(0,0) / min(o^[(l + l)e]^o(l,2/ / )^)^ , (3.9a) 

^( e )^( TO + 1 )" 1 ( 1 -^)^°'°) / Twn^lil-^SoMr^dj/ (3.96) 
These can be written as 



<(m + l)- 1 (l + h*? 4(0,0) [ min( - ' S (l,yT^dy' (3.10a) 
e d M JJ % [(l + jj)e]d 



^ > (m+ 1)-\1 - h^mO) [ min( 6 , ,S (1, y')^) dy' (3.106) 

N J J, 2[(l-^)e] d 



We now take limits of (3.10a) — (3.106) as e — > 0. We obtain 



li min f^!l>( m + l)-i(l_ 1)^0(0,0) / So(WrW (3-lla) 
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l imsup ^4-^ < (m + l)- 1 (l + ^-)^ i 0(0,0) / S Q (\,y')-* dy' (3-116) 
e^o ed JM J Ji 

Note that the integrals here are automatically finite since we are in superadapted coordi- 
nates and therefore all zeroes of So(l,y') are of order at most d — 1. We now take limits 
as N — > oo. If the endpoints of Ji are two zeroes z and z' of So(l,y'), then the interval 
will converge to [z,z']. Otherwise, the left endpoint of Ji may converge to zero and the 
right endpoint may go off to oo. In any event, the Ji goes to some (possibly unbounded) 
interval K i: and (3.11a) — (3.116) both converge to 

(m + l)"V(0,0) / S (l,yT kd dy' (3.12) 

The integral in (3.12) is finite since So(l,y') must have degree greater than d and has no 
zeroes of order d or greater since the coordinate system being used is superadapted. 

Furthermore, as N — > oo, the upper bounds of Lemma 2.1 and (2.27a) for Ig^e) 

and Is%{ e ) g° to zero. There is no issue of the constants appearing depending on iV due 
to the cutoffs' dependence on iV; the bounds of Lemma 2.1 and (2.27a) depend on \ \4>\\oo 

1+ (e) 

and no other properties of <fi. Hence we conclude that lim e ^o — ^ — is given by the sum 

of (3.12) over all i, along with their analogues in the other three quadrants. This gives 
exactly (1.6) (recall the terms where S(x,x m y') is negative gives no contribution) and we 
are done. 

Lastly, we prove (1.13a) of Theorem 1.3a. Either So(l,y) is positive for some 
y > or — So(l,y) is positive for some y > 0. Since the result we are trying to prove 
is symmetric in S and —S, without loss of generality we may assume that So(l,y) has 
this property. Then (1.6) says that Is,<j>{z) > As^e^ for sufficiently small e. This implies 
(1.13a) since I\g\ !( j,(e) > Is,</>(^) an d we are done. 

4. Case 2 proofs. 

We now assume S(x, y) is a smooth phase function in case 2 of superadapted 
coordinates with d > 1. Hence the bisectrix intersects N(S) at (d, d). As in the case 1 
proofs, we consider 7^(e) and J<t ^(A) and we will do some subdivisions of the domains 
of these integrals to prove the estimates and formulas. We proceed as follows. If (d, d) is 
the lower vertex of a compact edge e2, denote its equation by t\ + 7712^2 = «2- Where TV 
is a large natural number, fixed for now, we let A 2 = {(x, y) G [0, 1] x [0, 1] : x < -^y"^} 

and B 2 = {(x,y) G [0,1] x [0,1] : < x < Ny^}. Define C 2 = A 2 U B 2 . If (d,d) 

is not the lower vertex of a compact edge (i.e. (d, d) is on the vertical ray), then define 
C 2 = {(x,y) G [0, 1] x [0, 1] : x < y L }- Here L is large number to be determined by our 
future arguments. Similarly, if there is a compact edge e\ whose upper vertex is (d, d), we 
write its equation as t\ +mit 2 = ol\. We then let A\ = {(x, y) G [0, 1] x [0, 1] : y < jjX mi } 
and B x = {(x, y) G [0, 1] x [0, 1] : ±x mi < y < Nx m ^}. We then define C 1 = A 1 UB 1 . If 
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(d, d) is not the upper vertex of a compact edge, then define C\ = {(x, y) E [0, 1] x [0, 1] : 
y < x L }. 

In all cases, define D = [0, 1] x [0, 1] — (A 2 U B 2 ) We will see that the terms 
Is,4>( € ) an d Js,4>W dominate; the contributions from C\ and C 2 to the main term of the 
asymptotics can be made arbitrarily small as L — > oo. In fact, for the case when d comes 
from a compact edge, Lemmas 2.1 and 2.2 give that 

< CN ~ v " kd i J &( e )i < c ^ ( 4 - la ) 

I^WI < CN-*\\\-* |J* (A)| < C N \\\-$ (4.1b) 

Lemma 4.1. For sufficiently small e we have 

\IF,M <C|ln(e)|e* | J^(A)| < Cln |A||A|-4 (4.2a) 

Furthermore, for i = 1, 2 if Cj derives from the horizontal or vertical ray we have 

&)\ < j^jlH^ \J^(X)\ < ^-lnlAHAl^ (4.26) 

Proof: We start with the proof of (4.2a). We divide D = D\ U D 2 , where Di = {(x, y) G 
D : y < x m }. Here m is chosen such that if e\ exists, then m < mi, and if e 2 exists 
then m > m 2 . The estimates for Di and D 2 are proven the same way, so we restrict 
our attention to proving the estimates for Di. Taylor expand S(x,y) arount the origin as 
S(x,y) = Y. a <M, b<M s abX a y b + E M {x,y), where like in (2.19) for < a, (3 < M we have 

\d^E M (x,y)\ < C(\x\ M -« + \y\ M -?) (4.3) 

Correspondingly, we have 

d d y S{x,y)= Y, s' ab x a y b + d d y E M {x,y) (4.4) 

a<M, b<M-d 

The Newton polygon of dyS(x, y) has a vertex at (d, 0) and thus s' d0 ^ 0. If (d, d) is the 
lower vertex of some compact edge e 2 , then there is a compact edge e' of the Newton 
polygon of dyS(x,y) containing (d, 0) with equation ti + m 2 t 2 = d. Hence every nonzero 
s' ab x a y b appearing in (4.4) satisfies a + m 2 b > d and we can rewrite (4.4) as 

d d y s(x, y ) = s' do x d + s '^ a v b 

a<d, 0<b<M-d, a+m 2 b>d 

+ J2 s' ab x a y b + d d E M (x,y) (4.5) 

d<a<M, 0<b<M-d, (a,b)^(d,0) 
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If (d, d) is not the lower vertex of such a compact edge, (4.5) is still valid if we take the first 
sum to be empty. If the first sum is not empty, then since y < N~ Tn2 x m2 for all (x, y) G D, 
if (x, y) is in D then each term s' ah x a y h in the first sum of (4.5) is bounded in absolute 
value by N~ rn2b \s' ah \x a+m2b < N~ rn2b \s' ah \x d . Thus if iV were chosen sufficiently large, the 
absolute value of the whole first sum is less than ||s^ |a; d . 

Next, note that the absolute value of a given term s' ab x a y h in the second sum is at 
most \s' ab \x d (x + y). As a result, if the support of 4>(x, y) is sufficiently small, then for (x, y) 
in this support the absolute value of the second sum is also at most ||s^ |a; . Similarly, if 
the support of y) is sufficiently small, then by (4.2) and the fact that y < x m , for (x, y) 
in this support \d d E M (x, y)\ can also be assumed to be at most ll^ol 2 ^- Consequently, 
for such (x,y) in the support of 4> we can assume 

\d d S(x,y)\ >l\s' d0 \x d (4.6) 

Denote the vertical cross section of D at x by D x . By (4.6) and the measure version of 
Van der Corput's lemma, for each x in this range we have 

\{yeD x :0<S(x,y)<e}\<Ce*- 

x 

Also, \{y G D x : < S(x, y) < e}| is at most \D X \, which is at most x m if the support of 
is sufficiently small, which we may assume. Thus we have 

\{y G D x : < S(x,y) < e}\ < Cmin(x m ,e^-) 

x 

Consequently, integrating with respect to y first one has 

\I^)\< C f™< xm ^-) dx ( 4 - 7 ) 
Jo x 

i , i 
For small enough e, the quantities x m and €<* - are equal at xq = e (m+1 '> d , with x m smaller 

on the left and smaller on the right. Doing a computation gives 

rl ^ j 1 l ! i 

/ min(x m , e^-) dx = + — |ln(e)|e^ (4.8) 

J x' m + 1 (m + l)d l Wl v ; 

Hence \Is,j,( € )\ < C|ln(e)|e^ as desired. As for the Js^(X), by (4.6) the traditional Van 
der Corput lemma in the y direction gives 



Consequently, we have 



\[ e lXS{x ' v) (t){x } y)dy\<C\\\-^x- 1 (4.9) 
Jd x 

r 1 i 

\J^(X)\<C mm(N- m x m ,\\\-*-)dx (4.10) 
Jo x 
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This is exactly (4.8) with e replaced by |A| 1 . Thus instead of (4.8) for large |A| we get 
the estimate 

|J&(A)|<Cln|A||A|-3 (4.11) 

This completes the proof of (4.2a). Equation (4.26) is done the same way; the only dif- 
ference is that x m is replaced by x L . Equation (4.8) and its oscillatory integral analogue 
then give (4.26) and we are done. 

We now have proven the upper bounds for the smooth case: 
Lemma 4.2. Equations (1.12) and (1.156) hold. 
Proof. Add (4.2a) to (4.26) or (4.1a) - (4.16). 
Our next result is an analogue of Lemma 3.2. 

Lemma 4.3. There is a neighborhood U of the origin such that if <f)(x, y) is supported in 
U and 4>{x,y) is zero in a neighborhood of the origin, then 

|J 5l *(e)| < Ce^ \JsA^\ < C\\\~* (4.12) 

Proof. Fix some N. By (4.1a) - (4.16), the upper bounds of (4.12) hold for J^(e), 
^s^(^)' ^( e )' anc ^ ^s^W- Thus it suffices to prove these upper bounds for Ig^,(e) 
and Jg^ie), as well as -ff^(A) and «7^,(A) if they derive from the vertical or horizontal 
ray. These are all done basically the same way, so we restrict our attention to Ig^e) and 
J|^(e). As in the proof of Lemma 4.1, we divide D = Di U D 2 along the curve y = x m . 
The two pieces are done similarly, so we will only consider D\, the part where y < x m . 

Each vertical cross section {D\) x of the set D\ is the a subset of the interval 
[0,x m ]. Hence there is some 5 > such that on D, (j)(x,y) = for x < 5. Doing the y 
integration first we have 

I^WI <C J \{ye(D 1 ) x :0<S(x,y)<e}\dx (4.13) 

By (4.6), \dyS(x, y) \ is bounded below on x > 5. Hence by the Van der Corput lemma in the 

y direction, we have \{y € (Di) x : < S(x,y) < e}\ < Ce^ uniformly in x > 5. Inserting 
this back into (4.13) gives the desired bounds. For the oscillatory integral, one similarly uses 
the Van der Corput lemma in the y direction to get | J^ Dl ^ e lXS ^ x ' y ^ (j)(x,y) dy\ < C|A|~d 
uniformly in x > 5. Thus 

\JssW\ = I / ( / e lXS ^<P(x,y)dy)dx\ < C\\\~ k * (4.14) 

Js y (£>!)„, 
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These are the sought-after bounds for d|^(A) and we are done. 

We now proceed to the proof of the explicit formula (1.7). The general method- 
ology is similar to that of the case 1 arguments of section 3. If one writes = 0i + 02, 
where 0i is supported in a smaller neighborhood of the origin and 02 is zero on a neigh- 
borhood of the origin then by Lemma 4.3, limsup e ^ s '^ < ' t = limsup e ^ s '' Pl V and 

ln(e)e<2 ln(e)e<2 

liminf e ^o t = liminf e ^o s '^ 1 -r ■ So when proving (1.7) one can always replace by 

ln(e)ed ln(e)ed 

0! at will, regardless of how small the support of 0i is. 

We will show if one first chooses the parameter L of (4.26) sufficiently large, and 
then chooses 0i to be supported in a sufficiently small neighborhood of the origin, then 
the above limsup and liminf, added to their analogues from the other three quadrants, 
can both be made as arbitrarily close to the limit given in (1.7). To do this, in view 
of (4.26), it suffices to show that if the support of 0i is sufficiently small, the quantities 

jD jD 

lim sup e ^ s ' 4 ' 1 1 and lim inf e ^ s '% , can be made arbitrarily close to the appropriate 

ln(e)e'3 ln(e)e^ 

expression. For a fixed L we will find lower bounds for the limsup and upper bounds for 
the liminf. In doing so, we will choose the parameter M of the Taylor expansions in terms 
of L, and then the parameter iV of (4.1a) — (4.16) in terms of L and M. Analogous to 
in section 3 taking limits as L goes to infinity, both expressions will converge to the same 
limit. Adding this limit to its analogues in the other 3 quadrants will give (1.7). 

We start with the following lemma. As before, SddX d y d denotes the (d, d) term of 
the Taylor expansion S{x, y) = J2 a b s a bX a y h of S at the origin. 

Lemma 4.4. There are constants f3 and C depending on S{x,y), and a neighborhood U 
of the origin depending on S(x,y) and L, such that \S(x,y) — SddX d y d \ < \CL~ /3 x d y d \ on 

Dnu. 

Proof. Analogous to (4.5), we may Taylor expand 

S{x,y)-s d dx d y d = £ s' ab x a y b 

M>a>d, M>b>d, (a,b)^(d,d) 

+ s ab x a y b + s ab x a y b + E M (x,y) (4.15) 

a<d, M>b>d, a+m2b>ce2 M>a>d, b<d, a-\-mib>cti 

Here the second term is nonempty only if (d, d) is the lower vertex of a compact edge, and 
the third term is nonempty only if (d, d) is the upper vertex of a compact edge. 

The first sum can be made less than j^x d y d in absolute value by making the 
radius of U sufficiently small depending on M and S(x, y). If the second sum is nonempty, 
then the domain D is a subset of {(x, y) G [0, 1] x [0, 1] : < y < N~ m2 x m2 }. As a result 
if one changes coordinates from (x,y) to (x,y'), where y' = x m2 y, D becomes a subset 
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of D' = {(x,y) G [0, 1] x [0, 1] : < y < N~ m2 }. Observe that a given term s ab x a y b of 
the second sum becomes s a bX a+m2b (y') b . Since a + m^h > «2 and b > d in each term 
in the second sum, the entire sum can be written as y'(x a2 (y') d )f(x,y') for some f(x,y') 
which is a polynomial in y and a fractional power of x. Thus the sum is of absolute 
value at most CMN~ m2 x a2 {y') d in a small enough neighborhood of the origin. Note that 
Sdd,x d y d = SddX d+dm2 (y') d , and this is equal to SddX a2 (y') d since (d, d) is on the edge with 
equation t\ + m^ = 0.2- As a result, in the original (x, y) coordinates, the sum is of 
absolute value at most CMN~ m2 x d y d . Thus if one chooses N sufficiently large for fixed L 
and M, one has the desired bounds. 

The third sum is dealt with in exactly the same way, reversing the roles of the x 
and y axes. Since D necessarily lies in the range xl > y > x L , the error term EM(x,y) 
can be made less than j^x d y d by making the radius of U sufficiently small. This completes 
the proof of Lemma 4.4. 

Proof of (1.7). As before it suffices to assume 0(0,0) > as the general case can be 
obtained by writing = 0' — 0" where 0'(O,O) and 0"(O,O) > are positive. Let 8l > 
be such that on the ball -6(0, Sl) one has 

(1 - L-P)<f>(0, 0) < 0(x, y)<(l + L-P)<f>(0, 0) (4.16) 

Further assume that Sl is small enough that B(0, 8l) C U, where U is as in the previous 
lemma. Let ip(x,y) be a nonnegative cutoff function such that < ip < 1, tp(x,y) is 
supported on B(0, Sl), and ip(x, y) = 1 on [0, 4^] x [0, 4^]. Then by the discussion following 
(4.14) we may replace <f>(x, y) by <j)(x, y)ip(x, y) without affecting the liminf or limsup. We 
have 

I^(e) < (1 + L-^)0(O,O)|{(x,y) G D : < S(x,y) < e}\ 
By Lemma 4.4 this is bounded by 

(1 + L"^)0(O, 0) \{(x, y) e D : < (s dd - CL^)x d y d < e}\ (4.17a) 

In addition, 

isA*) X 1 " L ~ P )^ °M X > v) e D n y] x [°» y]) : < s ( x > v) < e >l 

> (1 - L"' 3 )0(O, 0)|{(x, y) E D n [0, (^] x [0, S f}) : < (s dd + CL^)x d y d < e}\ (4.176) 

If s dd < 0, then by (4.17a) Ig ^(e) = for large enough L. For s dd > 0, we need the 
following lemma, whose proof is routine. 

Lemma 4.5. Suppose < m 2 < mi and <5 > 0. Then as t — > 0, 

|{(x,y) G (0,5 ] x (0,6o] :x m ^ <y<x m2 , y< -}\ = (-L- ^-)tln(t) + 0(t) 

1 x mi + 1 mo + 1 
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|{(x, y) E (0, So] x (0, So] : < y < x m \ y < -}| = ^-j-tln(t) + 0(t) 

We now apply Lemma 4.5 to (4.17a) - (4.176). We get that 

Jf» < (1 + L-P)(s dd - CL-Py^iO, 0)e^ ln(e)(^— L-) + Ce^ (4.18a) 

^ mi + 1 m 2 + 1 

> (1 - L^)(s dd + CL-f>)-i<K0, 0)6^ ln(e)(-^— - ) - Ce^ (4.186) 

(When (d, d) is on the horizontal ray one substitutes m\ = L, and when it is on the vertical 
ray one substitutes m 2 = 7;)- Hence we have 

hmsup < (1 + L-P)(s dd - CL-*)-40(O, 0)(^— L-) (4.19a) 

e^o ed ln(e) mi + 1 m 2 + 1 

liminf > (1 - L^)( Sdd + CL^)-*0(O, 0)(^— - -J— ) (4.196) 

C d ln(e) mi + 1 m 2 + 1 

We take limits as L — > 00. Both expressions converge to 0(0, 0)( mi 1 +1 — m 1 +1 ) ) where now 
mi is taken as 00 when (d, d) is on the horizontal ray and m 2 is taken as when it is on 
the vertical ray. Hence by the discussion following (4.14) we conclude that 

itAe) 1 1 1 
lim ^*1L = S 7/0(O,O)( -) 

Letting Si be the slope — —of the edge t\ + mfo = oti, this becomes 

Till 



itAe) 1 1 1 



lim f'^ v = s 
e~d ln(e) 



dd 



0(0, 0)(- -) (4.20) 

Si — 1 S2 — 1 



In summary, if s dd > then (4.20) gives the contribution to (1.7) from the upper right- 
hand quadrant. If s dd < then the contribution is zero as mentioned above Lemma 4.5. 
Adding this to its analogues over the other three quadrants gives exactly the formula of 
(1.7) and we are done. 

Our final task is to prove (1.136): 

Proof of (1.13b). Since the result is symmetric in S and — S, we may replace S by —S 
if necessary and assume that s dd > 0. As in the proof of (1.13a), we write = 0i + 2 , 
where 0i is nonnegative and 2 is zero on a neighborhood of the origin. By Lemma 4.3, 

lim e ^ J|s| ^ 2 i' > = 0. So to prove (1.136) it suffices to show Ji^i < ^ 1 (e) > As^il ln(e)|e^ for 

ln(e)e 3 ' 

some As, fa > 0. Since I\s\,fa(£) > Is,fa{t)i it further suffices to show that Is, fait) > 
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As^l ln(e)|e^ for some As, fa. For this we use Lemma 4.4, which implies that there is a 
5 > such that on ([0, 5] x [0,5]) n D we have 

S{x,y) <^s dd x d y d (4.21) 

As a result, shrinking 5 to ensure that <j)(x,y) > |0(O,O) on [0,5] x [0,(5] if necessary, we 
have 

W<0 > |0(O, 0)K(x,y) 6DH ([0,5] x [0,5]) : ^ dd xV < e}| (4.22) 

Using Lemma 4.5, we conclude that there is some As, fa with (is,^ (x, y)| > Ag^J ln(e)|e^ 
as needed. This gives (1.136) and we are done. 

5. Case 3 proofs. 

In this section, S(x,y) is a smooth phase function in case 3 of superadapted 
coordinates with d > 1. We restrict ourselves to the situation where the bisectrix intersects 
the horizontal ray in its interior, as the case of a vertical ray is entirely analogous. Thus 
the lowest vertex of N(S) is of the form (c,d), where c < d and d is also the Newton 
distance of S. As in sections 3 and 4 we will focus our attention on the analysis of Is<t>( 6 ) 
and Jg^iX). We divide [0, 1] x [0, 1] into two parts. For a sufficiently large positive integer 
k (to be determined by our arguments), we let D\ = {(x, y) E [0, 1] x [0, 1] : y < x k } and 
D 2 = {(x, y) E [0, 1] x [0, 1] : y > x k }. Our first lemma is the following. 

Lemma 5.1. |if^(e)| < Ce^ and |j£j(A)| < C\X\~^. 

Proof. As in cases 1 and 2, we write the Taylor expansion of S at the origin as S(x, y) = 
^Za<M b<M s abX a y h + EM(x,y), where for < a, (3 < M the function Em(x,u) satisfies 
the error estimates 

\d^ y E M {x,y)\ < C(\x\ M - a + \y\ M -?) (5.1) 
The dth y-derivative can be written as 

d d S(x,y) = s' ab x a y b + d d E M (x,y) (5.2) 

a<M, b<M-d 

Furthermore the Newton polygon of d d S(x,y) has a vertex at (c, 0), contained either in 
the vertical ray of N(d d S) or an edge of N(d d S) with equation t\ + mti = c with m > 0. 
Hence each (a, b) in the sum of (5.2) satisfies a + mb > c. Analogous to (4.5) we rewrite 
(5.2) as 

d d S(x,y) = s' c0 x c + <bX a y b 

a<c, 0<b<M — d, a+mb>c 

+ J2 s' ab x a y b + d d E M (x,y) (5.3) 

c<a<M, 0<b<M-d, (o,6)^(c,0) 
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In the case where (c, 0) is on a vertical ray of N(dyS) : the first sum of (5.3) is empty. 
We now argue like after (4.5). Since y < x k for all (x,y) G Di, if (x,y) is in Di 
each term s' ah x a y h in the first sum of (4.5) is bounded in absolute value by \s' ab \x a+kb < 
\s ab \x k ~ m (x a+mb ) < \s' ab \x k ~ rn+c . Thus as long as k were chosen greater than m, which 
we may assume, then if the support of y) is sufficiently small, then for (x, y) in this 
support, the absolute value of the whole first sum is less than ||s' c0 |^ c - Also, the absolute 
value of a given term s' ab x a y b in the second sum is at most \s' ab \x c (x + y). As a result, for 
such (x,y) the absolute value of the second sum is also at most lls^o! 3 ^- Similarly, using 
(5.1) and the fact that < y < x k , for such (x,y) the quantity \dyEM(x,y)\ can also be 
assumed to be at most ||Scol xC - Consequently, for these (x,y) we have 

\d d y S{x,y)\>- A \s' c0 \x c (5.4) 

As a result, by the measure version of Van der Corput's lemma of [C], for each x in this 
range we have 

\{y : < S(x, y) < e}| < Ce^x~^ (5.5) 
Consequently, integrating with respect to y first one has 

\Isd>( e )\< C f e*x-% dx = C'e* (5.6) 
Jo 

This is the desired upper bound for Is\{ e )- As for the oscillatory integral analogue, by 
(5.4) the normal Van der Corput lemma in the y direction gives 

Hence by integrating first with respect to y one has 

\J&(\)\<C [ 1 \\\-*x-*dx = C , \\\-* (5.8) 
Jo 

This completes the proof of Lemma 5.1. 

It turns out that one gets stronger estimates for the Ig 2 ^(e) and jf^A). Observe that 
since c < d, the quantity is greater than |. We have the following. 

Lemma 5.2. |/^(e)| < Ce^ and |J^(A)| < C\X\~^. 

Proof. We will verify the hypotheses of Lemma 2.0, with the roles of the x and the y 
variables reversed. Because (c,d) is the rightmost vertex of N(S), the Newton polygon of 
d°S(x,y) has a single vertex at (0, d). The Taylor expansion of d°S(x,y) can be written 
in the form 

d c x S(x, y) = r od y d + J] r ab x a y b + E M (x, y) (5.9) 

0<a<M, d<b<M (a, 6)^(0, d) 
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Similar to elsewhere in this paper, bounding the error term using the fact that x < y * on 
D 2 , in a small enough neighborhood of the origin on D 2 one has 

\d c x S{x,y)\> 1 -\r Qd \y d 

Thus if c > 2, one can apply Lemma 2.0 and immediately get this lemma. If c = 1, to 
apply Lemma 2.0 one also needs that \d^S(x,y)\ < Cy d ~k. But in fact since the Newton 
polygon of d%S(x,y) is a subset of {(#, y) : y > d} : by expanding as in (5.9) one even 
has the stronger estimate \d%S(x,y)\ < Cy d . Thus Lemma 2.0 applies here. If c = 0, to 
apply Lemma 2.0 one needs (2.2c) to hold (with the x and y variables reversed) which 
here means one needs \d y S(x,y)\ > Cy d ~ x and \dyS(x,y)\ < C'y d ~ 2 . Since the Newton 
polygon of dyS(x, y) has a single vertex at (0, d — i) this holds as in the c = 1 case. Lastly, 
to apply Lemma 2.0 for c = one also needs that d > | + 1. We can make this true simply 
choosing k large enough since d is at least 2. This completes the proof of Lemma 5.2. 

Corollary 5.3. (1.11) and (1.15a) hold in Case 3. 

Proof. Add the estimates from Lemmas 5.1 and 5.2 and their analogues from the other 
three quadrants. 

Next, we prove the lower bounds of (1.13c). Since we are not trying to prove 
sharp estimates, the arguments are not that intricate. Assume 0(0, 0) 7^ 0, and let M 
be some large positive integer. We examine the behavior of S(x, y) on the set Zn = 
{(x,y) : x > 0, x N < y < 2x N } for iV sufficiently large. We Taylor expand S(x,y) = 
12 a <M b<M s abX a y h -\-Em{x, y) as above. If N is large enough, the term s c dX c y d dominates 
this Taylor expansion much the way r$dy d dominates (5.9) or s' c0 x c dominates (5.3). Hence 
in a small enough neighborhood U of the origin, on we have 

\S(x,y)\<2\s cd \x c y d (5.10) 
Shrinking U, we may assume that \(p(x,y)\ > ±|<£(0, 0)|. Hence for C<t> = ±|0(O,O)| we have 
I\s\A e ) > c ^l{( x ' y)£UnZ N : 2\s cd \x c y d < e}\ (5.11) 

It is easy to compute that the curve y = 2x N intersects the curve 2\s c d\x c y d = e at 
x = ae^fd+c f or some a depending on S(x,y). Hence for e sufficiently small the measure of 
the set in the right-hand side of (5.11) is at least the measure of the portion of Zn between 

1 1 N + l 

x = |e JVd + c and x = ae Nd + c , given by a e Nd + c where now a also depends on N. Thus we 
can write 

I\s\A e ) > c We^ (5-12) 

Since d > c, the exponent in (5.12) is larger than ± but as N — > 00 it tends to |. This 
gives us (1.13c). 
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We now move to the case of real-analytic phase. Our goal here is to prove Theorem 
1.1 c). So assume S(x,y) is real-analytic. It suffices to show that lim e ^ Js '^ exists and 



T+ 



1 (e) 

is given by (1.8). As in cases 1 and 2, we will give an expression for lim e ^ -^4 — and 
the full limit will follow by adding this and the analogues from the other quadrants. Also, 

I° 2 (e) i Dl (e) 

by Lemma 5.2, lim e ^o ^4 — = so it suffices to show lim e ^o -^4 — exists and has the 

desired value. Because the bisectrix intersects N(S) in the interior of its horizontal ray 
and (c, d) is the lowest vertex of N(S), the real-analytic S(x, y) can be written as 

S(x,y) = s cd x c y d + x c+1 y d g(x) + S «* X V (5.13a) 

b>d+i 

Here g(x) is real-analytic. Changing coordinates from (x, y) to (x, y') where y = x k y', we 
have 

S(x,x k y')=s cd x c + kd (y') d + x c + kd + 1 (y') d g(x)+ £ s ab x a+kb (y') b (5.136) 

b>d+i 

Since the line t± + rati = a is an edge of N(S) containing (c, d), each (a,b) in the sum 
(5.13a) satisfies a + mb > c + md. Furthermore, b > d, and therefore 

a + kb = a + mb + (k — m)b > c + md + (k — m)b > c + md + (k — m)d = c + kd 

Hence we can rewrite J2b>d+i s abX a+hb (y') b as x c+kd+1 (y') d+1 f(x, y) where f(x, y) is real- 
analytic. Thus we have 

S(x,x k y') = s cd x c+kd (y') d + x c+kd+1 (y') d g(x) + x c+kd+1 (y') d+1 f(x,y') (5.14) 

In the (x,y') coordinates Is\{ e ) becomes 

I^Je)= [ x k cf>(x,x k y')dxdy' (5.15) 

i{(x,l(')e[0,l]x[0,l]:0<S(i8,xV)<e} 

Let d! be between c and d. The exact value of d! will be dictated by our arguments. Then 

i 

1 £ d'(fc + l) , J_ 

the portion of (5.15) over x < e d '( fe+1 ) has absolute value at most C J Q x dx = Ce> 



Since this is o(ej), this portion of the integral will can be removed without affecting 

I Dl (e) 

lim e ^o -^4 — • In other words, we may replace -^^(e) by /^(e) where 



I's,4e)= [ , x k cf>(x,x k y')dxdy' (5.16) 

J {{x,y')e[e d '( k + 1 ) ,l]x[0,l]:0<S(x,x k y')<e} 

We now fix x > t d ' ( - k + 1 ) and look at the set E x = {y e [0, 1] : < S(x,x k y') < e}. We 
may assume the support of (f)(x,y) is small enough so that if 4>{x,y) ^ then x is small 
enough so that 

\x c+kd+1 (y') d g(x)\ + \x c+kd+1 (y') d+1 f(x,y)\ < \\s cd \x c+kd y d (5.17) 
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Note that if s c d is negative, by (5.14) and (5.17) S(x,x y') is always negative and thus 
I' s ^(e) = 0. So assume that s c a > 0. Then by (5.14) and (5.17) we have 



E x C {y G [0, 1] : < s cd x c+fc V < 2e} C [0, Ce^x"^] (5.18) 
i 

Since x > e d '( fe+1 ) , we have 

1 c+fcd 1 c + fcd fcd' + d' c + fcd k(d' ~d) + (d' -c) 

a - < e d 6 d(d')(fe+i) = e d(d')(fe+i) e d(d')(fe+i) = e d^fe+T) 

i 

Hence if d! were chosen close enough to d, there is some r)' > such that for x > e d '( fe+1 > 
one has e^x ~ < e v and thus 

£ x C [0,Ce"'] (5.19) 
Next, we write J^(e) = 7^(e) + /^(e), where 

7^(e) = / i x k (f)(x,0)dxdy' (5.20a) 

I'^(e)= f , x fc (0(x,xV)- 0(x,O))dxdj/' (5.206) 

^{(^,y')e[e d ' (fe + 1) ,l]x[0,l]:0<S(x,x fe y')< £ } 

Note that due to (5.19), the factor (<f)(x,x k y') — 0(x,O)) in (5.206) is bounded in absolute 
value by C'e 71 , so we have 

I's,4e)<C'e v ' [ x k dxdy' (5.21) 

In (5.21), we perform the y integration by inserting the second inclusion of (5.18). We 
then have 

I'sAe) < C"e"' ^ , x k e^x~ £± ^ dxdy' < C"'e^' /V§ dx = C""e* +r >' (5.22) 
J £ d'(k+i) J 

Thus lim e ^o -~r — = 0. Hence lim e ^o _J ^4 — = lim e ^o -~r — ; and our goal now becomes 

£ d £ d £ d 

to prove the latter limit gives the portion of (1.8) coming from the upper right quadrant. 
Next, we rewrite 1$ ^(e) as 

I^(e) = f_ , x k ( f>(x,0)\{y' : < S(x,x k y') < e}\ dxdy' (5.23) 

J e d'(k+i) 

To analyze (5.23), first note that - dX c+kd^ x c+kd+i g ^ is a real-analytic function in a neigh- 
borhood of the origin, which we denote by h(x). Then by (5.14) we have 

S(x,x k y') = (s cd x c+kd + x c+kd+1 g(x))((y') d + xh(x)(y') d+1 f(x,y')) (5.24) 
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Since s cd is being assumed to be positive, in (5.24) (s cc ix c+kd + x c+kd+1 g(x))^ is positive 
and there is j(x, y') with d y j(x, 0) = 1 and d x j(x, 0) = such that (5.24) can be rewritten 

as 

S(x,x k y')$ = (s cd x c+kd + x c + kd+1 g(x)^j(x,y') (5.25) 
Consequently, in (5.23), one has 

\{y' : < S(x,x k y') < e}\ = \{y' : < j(x,y') < e^{s cd x c+kd + x c+kd+1 g(x))-^}\ (5.26) 

By the inverse function theorem, (x,j(x, y')) has an inverse function which can be written 
as (x, k(x, y')) for some k(x,y') which satisfies d y k(x,0) = 1 and d x k(x,0) = By (5.19) 
the interval of (5.26) has length at most Ce v . As a result, as long as e is small enough we 
can use a linear approximation to k(x, y') and get that 

\{y' : < j(x,y') < e^{s cd x c+kd + x c+kd+1 g(x))-^}\ 

= (s cd x c+kd + x c+kd+1 g(x))-^ + e"'0(e* (s cd x c+kd + x c+kd+1 g(x))~*) (5.27) 
Thus we have 

I^(e) = I , x k <P(x,0)eHs cd x c+kd + x c+kd+1 g(x))-Ux 

+0{e r >' / i x k \<j)(x,0)\e^(s cd x c+kd + x c+kd+1 g(x))~^ dx) (5.28) 

Because rj' > 0, the ratio of the second term to goes to zero as e > (assuming the 
integral is finite, which we will see shortly). Thus the second term does not contribute to 

lim e ^ s ' 4 'i 6 and we have 

I" (e) f 1 

lim^4^= / x k (s cd x c+kd + x c+kd+1 g(x))- 1 d ( j)(x,0)dx 
ed Jo 

= [ (s cd x c + x c+1 g(x))-^(p(x,0) dx (5.29) 
Jo 

Since \s cd x c+kd + x c+kd+1 g(x) \ > ^x c+kd and c < d the integral (5.29) is finite as needed. 
Going back to the definition of g{x), S{x,y) = (s cd x c + x c+1 g(x))y d + 0(y d+1 ). As a 
result, (5.29) translates into the part of equation (1.8) coming from the upper right-hand 
quadrant. (Recall that (5.29) is for s cd > and that the limit is zero when s cd < 0). 
Adding the analogous expressions from the remaining three quadrants gives (1.8) and we 
are done. 

6. Proof of Theorem 1.6b. 
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Suppose S(x,y) is a smooth phase function in case 1 superadapted coordinates, 
and 4>(x, y) is nonnegative with 0(0,0) > 0. Let ijj(t) be a nonnegative function in C C (R) 
such that ip(t) > 1 on [—1, 1]. Then by (1.13a), if j is sufficiently large we have 

| / i/>(2iS(x,y))<l>(x t y)dxdy\ > | / (f>(x,y) dx dy\ > A s ^2~^ (6.1) 

J J{(x,y):\S(x,y)\<2-i} 

We also have 

' ^{2 J S{x, y))(p(x, y) dx dy = 2~ j I ( / ^ (2^ ' X)e lXS(x ' y) dX)dx dy 



= 2~ j J *p(2- J \)( J e lXS(x ' y) (j)(x,y)dxdy)dX = 2~ J J ^(2~ J A) JsA X ) dX ( 6 - 2 ) 

In order to prove (1.16a), we argue by contradiction. Suppose that we were in the setup 
of (1.16a) but limsup A ^ 00 I Js '*\ x > I = 0. Then for any 5 > 0, we may let Ms be such that 



A d 



for |A| > Ms we have \Js,<t>{X)\ < 5\\\ 3. We then have 

\2~i j ^(2^ A)J S)0 (A) dX\ = 



|2" J / $(2-i\)J s ,4,{\)d\ + 2-i [ V>(2^'A) J S , (A) dX\ 

J\\\<M S J\\\>M S 

< C2~iM 5 + 2 - J 5 [ |^(2-''A)||A|-3dA (6.3) 



In turn, equation (6.3) is bounded by 

C2- J Ms + 2~ J 5 [ |^(2- J A)||Ar^dA 
Jr 

= C2- j M s + 2-U I \i>{X)\\X\~^dX 
Jr 

<C2- J M 5 + C'52-i (6.4) 

When j is sufficiently large, (6.4) is at most 2C'52~i. On the other hand, by (6.1), it 
must also be at least As^2~i. This gives a contradiction if 5 were chosen less than 



This contradiction implies that limsup^^ | Js -^ A - ) | is in fact positive, giving (1.16a). 

A d 

That the limsup of (1.16c) is positive is proven from (1.13c) exactly as (1.16a) 
is proven from (1.13a), so we do not include a proof here. Since it holds for all S > 
the limsup is automatically infinite. Equation (1.166) is proved similarly to (1.16a), using 
(1.136) in place of (1.13a). Namely, suppose S(x,y) is a smooth phase function in case 2 
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superadapted coordinates, and (f>(x, y) is a nonnegative function with 0(0, 0) > 0. Using 
(1.136) we have 



| / ^(yS(x,y))<P(x,y)dxdy\ > | / <f>{x,y) dx dy\ > A s J-2~^ (6.5) 

J J{{x, V y.\S{x,y)\<2-i} a 

Exactly as above we also have 

J iP(VS(x,y))ct)(x,y)dxdy = 2-i J ^(2"'A) J^(A) dX (6.6) 

Proceeding by contradiction again, suppose (1.166) does not hold. Therefore for every 
5 > there is some L$ such that for |A| > L$ we have \ Js,4>(X)\ < S\X\~d In |A|. Analogous 
to (6.3) we have 

2~ j \ [ 4>(2- j X)J s ,4>( X ) d M < C2~ j L s + 2~ J 5 [ |^(2" J A)| \X\~hn \X\ dX 

< C2~ j L s + 2~ j 5 [ |V>(2- J A)||A|-^ln|A|(iA (6-7) 
Jr 

Changing variables, this in turn is equal to 

C2- j L s + 2~iS [ |^(A)||A|-* ln(2 J ' | A|)cZA 



'R 



= C2~ j L s + 2-U I \^{X)\\X\- k d\n\X\dX + j2-i\n{2)8 I |^(A)||A| _ a dA 
Jr Jr 

<C2- j L s + C"6j2-i (6.8) 

If j is sufficiently large, (6.4) is at most 2C"5j2~i , while by (6.5) it is at least As,<j,^2~i. 
This is a contradiction if 5 < 4^4t- Therefore limsup;,^^ I Js '^ x \ I must in fact be 

2dC l ln(A)A~d l 

positive. Hence we have (1.166) and we are done. 

7. Superadapted coordinates. 

In this section we prove the existence of superadapted coordinates for smooth 
phase functions. Here we always assume S(x, y) is a smooth phase function defined on 
a neighborhood of the origin such that 5(0,0) = and S(x,y) has nonvanishing Taylor 
expansion at the origin. The three cases of superadapted coordinates can be written as 
follows: 

Case 1. The bisectrix intersects N(S) in the interior of a bounded edge e and any real 
zero r ^ of S e (l,y) or S e (— 1, y) has order less than d(S). 
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Case 2. The bisectrix intersects N(S) at a vertex (d, d) and if e is a compact edge of 
N(S) containing (d, d) then any real zero r 7^ of S e (l, y) or S e (— 1, y) has order less than 
d(S). 

Case 3. The bisectrix intersects N(S) in the interior of one of the unbounded edges. 
Lemma 7.0. Any superadapted coordinate system is adapted. 

Proof. By the main theorem of [Gl], if U is a small enough neighborhood of the origin 
and eo denotes the supremum of the numbers e for which $ v \S\~ e is finite, then d(S) < — , 
with d(S) = ^ in cases 1, 2, and 3. Hence if one is in cases 1, 2, or 3, one is in adapted 
coordinates. 

Case 2 has some special features for which the following preliminary lemma will be useful. 
Related lemmas occur in [PSSt] and [V]. 

Lemma 7.1. Suppose the bisectrix intersects N(S) at a vertex (d,d) but is not in super- 
adapted coordinates. Correspondingly, let e be a compact edge of N(S) containing (d, d) 
such that S e (l, y) or S e (—1, y) has a zero of order d or greater. If (d, d) is the upper vertex 
of e, then S e (x, y) is of the form cx a (^- — r) d for positive integers a, m and some nonzero 
c, r. If (d, d) is the lower vertex of e, then S e (x, y) has the analogous form c'y a (-^r — r') d . 
with c', r' 7^ and a', m' positive integers. 

Proof. We first consider the case where (d, d) is the upper vertex of e. Write the equation 
of e as t\ +mt 2 = a. We will show that these values of m and a work. Note that if s a bX a y b 
appears in S e (x,y) then a + mb = a. We factor out x a , writing S e (x,y) = x a T e (x,y). 
Each term of T e (x, y) is now of the form t a bX a ~ OL y h with (a — a) +mb = or (a — a) = —mb. 
Thus we have 

t ab x a - a y b = U(4r) 6 ( 7 -!) 
Conequently for a polynomial P(z), we can write 

S e (x,y)=x a P(^) (7.2) 

Plugging in x = 1 or —1, we see that P(y) is a polynomial of degree d with a real zero r 
of order d or greater. Therefore we must have P(y) = c(y — r) d for some c ^ 0. Hence 
S e (x, y) = cx Q (^- — r) d . Since x can only appear to integer powers, m must be an integer 
and therefore a is as well. Also, since it comes from an edge S e (x,y) contains multiple 
terms. Hence r ^ and we are done with the case where (d, d) is the upper vertex of e. 

The case where (d, d) is the upper vertex of e is done similarly. Since S e (l,y) or 
S e (— l,y) has a zero r ^ of order d or more, S e (x,y) has zeroes of order d along some 
curve y = rx m . Hence S e (x, 1) or S e (—x, 1) has a zero not at the origin of order d or more 
and the above argument applies, reversing the roles of the x and y variables. 
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Lemma 7.2. Suppose one is not in superadapted coordinates. Suppose e is an edge of 
N(S) intersecting the bisectrix in its interior with equation t\ + rati = a for m > 1 such 
that S e (l,y) has a zero r 7^ of order k > d(S). Then m is an integer and both S e (l,y) 
and S e (— 1, y) have a zero of order k not at the origin. 

Proof. Exactly as (7.2), there is some polynomial Q(y) such that for x > we have 

Se (x,y)=x a Q(^-) (7.2') 
Plugging in x = 1, we see that Q(y) = S e (l,y). 

We now show that m must in fact be an integer. To see this, note that if m 
were not an integer, then the degrees of the powers of y appearing in S e (l,y) would have 
to be separated by at least 2. Hence S e (l,y) would have to be of the form y l3 R(y c ) for 
some (3 > 0, c > 2, where R is a polynomial. Next, since (d(S), d(S)) is on N(S), we have 
a = (1 + m)d(S). Since m > 1 when m > 1 is not an integer, the maximum possible value 
of yon the line t 1 + mt 2 = (l + m)d(S) for ti,* 2 > is m ± k d(S) < 2d(S). Thus the degree 
of y@R(y c ) is less than 2d(S), and hence the degree of R(y) is less than 2d ^ < c^S 1 ). Hence 
the zeroes of -R(y) are of order less than d(S) : implying the zeroes of of S e (l, y) = y f3 R(y c ) 
other than y — are of order less than d(S). This contradicts our assumption that S e (l,y) 
has a zero r 7^ of order k > d(S) and we conclude that m is an integer. 

Note that since m is an integer so is a. Consequently by (7.2') if m is even, then 
S e (l,y) = ±S e (— while if m is odd one has S e (l,y) = ±S e (—l,—y). Hence in either 
case both S e (l,y) and S e (— l,y) have a zero of order & not at the origin. This completes 
the proof of Lemma 7.2. 

The next lemma is the crux of this section. To set it up, suppose S(x, y) is not 
in superadapted coordinates and the bisectrix intersects the interior of an edge e. Then 
since «S e (l,y) or S e (—l,y) has a zero r 7^ of order k > d(S), S e (x,y) has zeroes of order 
k at any point on a curve of the form y = rx m . Hence S e (x, 1) or S e (x, —1) has a zero 
of order k away from the origin. Thus we may switch the roles of the x and y axes if we 
want and assume e has equation t\ + mti = a for m > 1; by Lemma 7.2 m is an integer 
and S e (l, y) has a zero r 7^ of order k. 

Lemma 7.3. Suppose S(x, y) is not in superadapted coordinates and the bisectrix in- 
tersects the interior of an edge e. As described above, switching the x and y axes if 
necessary, write the equation of e as t\ + mt 2 = a for an integer m > 1 and assume 
S e (l,y) has a zero r 7^ of order k > d(S). Then there is a coordinate change of the form 
(x,y) — > (x,y + a(x)) such that a(x) is smooth with a(0) = 0, after which one is either 
case 1 or 3 of superadapted coordinates, or the following more general version of case 2: 

Case 2'. The bisectrix intersects N(S) at a vertex (d,d). 

Proof. Let Q(y) = S e (l,y), and let (p, q) denote the upper vertex of the edge e; necessarily 
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q > d(S). We will find a smooth function a(x) such that S'(x, y) = S(x,y + a(x)) is in one 
of the following two mutually exclusive categories. 



Category 1: S'(x,y) is either in case 1, case 2', or case 3. 

Category 2: The bisectrix intersects the interior of an edge e' of N(S') with equation 
ti + m't2 = ol , m' > m > 1, such that the upper vertex (p',q') of e' satisfies q' < q and 
such that S'(x,y) is not in case 1. (In particular by Lemma 7.2 S' e , (l,y) has a zero of 
order > d(S')). 

Lemma 7.3 will then follow; there can be at most q iterations of category 2. 

We first consider the case where k < q. The function Q(y + r) has a root at y = 
of order k. We choose a(x) = rx m and define 5"(x, y) = S(x, y + a(x)) = S(x, y + rx m ). 
Note that ti + mt 2 = a is a supporting line of N(S') as it was for N(S), and that 
there is an edge E of N(S') on this line whose upper vertex is (p,q)- Observe that 
S' E (x, y) = S e (x, y + rx m ) = x a Q{-^ + r). Since Q has a zero of order k at r, the lowest 
power of y appearing in S' E (x,y) is y k and therefore _E"s lower vertex is at a point (j, k) 
for some j. Since both vertices of E have y-coordinates at least d(S), they are both in the 
portion of the line ti + m£ 2 = a on or above (d(S), d(S)). Thus the edge E lies wholly on 
or above the bisectrix. If the bisectrix intersects N(S') at a vertex or inside the horizontal 
or vertical rays, one is in Category 1. Otherwise, it must intersect N(S') in the interior of 
an edge e' whose upper vertex is either (j, k) or a lower vertex. And because t± + mt<i = a 
is a supporting line for N(S') and e' lies below E, e' will have equation ti + m't2 = ol 
for some m' > m > 1. Thus we are either in case 1 superadapted coordinates (which is in 
Category 1) or we are in Category 2. Hence when k < q, S'(x,y) is in either Category 1 
or 2 and we are done. 

It remains to consider the situation where r is a zero of Q(y) of order q. In this 
case we have Q(y) = c(y — r) q for some c. For a large integer n we expand S(x, y) as 

S(x, y) = cx a (^- rf + T n (x, y) + E n (x, y) (7.3) 

Here the polynomial T n (x,y) are the terms of S"s Taylor expansion with exponents less 
than n. For all < j3, 7 < n one has 

' n ix,y)\<C(\x\ n -e + \ y \ n -<) (7.4) 



dx^dy" 



Note that 



S(x, x m y) = cx a {y - r)« + x a+1 T^(x, y) + E n (x, x m y) (7.5) 
Here T^(x,y) is also a polynomial. We define s(x,y) = ^ ^ , so that 

s(x, y) = c(y - rY + y) + x~ a E n (x, x m y) (7.6) 
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We claim that the function s(x,y) is smooth on a neighborhood of (0,r). Off the y-axis 
smoothness holds because S(x, y) is smooth. One can show that a given derivative of 
s(x,y) exists when x = and equals that of c(y — r) q + xT^(x,y) for large enough n 
by examining the difference quotient of a one-lower order derivative of (7.6), inductively 
assuming this lower-order derivative exists and has the right value when x = 0. Equation 
(7.4) ensures that the difference quotient of the lower derivative of x~ a E n (x,x m y) tends 
to zero as x goes to zero. We conclude that s(x,y) is smooth on a neighborhood of (0,r). 

We next use the smooth implicit function theorem on g 9 yq -t and find a smooth 

function k(x) defined in a neighborhood of x = such that k(0) = r and q 9 q -t (x, k(x)) = 0. 
Transferring this back to S(x,y) we have 

d q ~ 1 S 

-^ ZT (x,x m k(x)) = (7.7) 

Thus if we let a(x) = x m k(x) and S'(x,y) = S(x,y + x m k{x)), for all x we consequently 
have 

d q ~ 1 S' 

^FT-(^ 0)=0 (7-8) 
Thus for every a the Taylor series coefficient S' aq _ 1 is zero. 

Next, since ti+mt 2 = a> is a supporting line for N(S), this line is also a supporting 
line for N(S') and intersects N(S') at the single vertex (p,q). If S'(x,y) is in Category 1 
we have nothing to prove, so we may assume we are not in Category 1. Let e! denote the 
edge of N(S') intersecting the bisectrix and denote its equation by t\ + m't^ = a' . Since 
e' lies within the set t\ + rati > ol and is no higher than the vertex (p, q) of N(S') that is 
on the supporting line t\ + rati = a, we have m' > m > 1. If the upper vertex (p' , q') of 
e' satisfies q' < q, one is in Category 2 and we are done. So we assume this upper vertex 
is (p, q) itself. 

If S' e , (l,y) has a real zero r' ^ of order k < q, one is in the situation above 
(7.3); there is a smooth b(x) such that S'(x, y + b(x)) = S(x,y + a(x) +b(x)) is in Category 
1 or 2 as needed. The only other possibility is that S' e ,(l,y) has a single zero r' ^ of 
order q. But this cannot happen. For this would imply S' e ,(x,y) = c'x a (-^r — r') q has 
a nonvanishing y q ~ x term. Consequently, for some a the Taylor series coefficient S' aq _ l 
would be nonzero, contradicting (7.8). Thus the case where S' e ,(l,y) has a single zero of 
order q does not occur, and we are done with the proof of Lemma 7.3. 

The final step of the proof of the existence of superadapted coordinates is the 
following. 

Lemma 7.4. Suppose one is in case 2'; that is, the bisectrix intersects N(S) at a vertex 
(d : d). Then there exists a smooth coordinate change fixing the origin after which one is 
in case 2 of superadapted coordinates. 
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Proof. Suppose the bisectrix intersects N(S) at a vertex (d, d) but S(x, y) is not in 
superadapted coordinates. Then (d,d) is on an edge e of N(S) such that S e (l,y) or 
S e (— 1,?/) has a zero r ^ of order o? or greater. By Lemma 7.1, S e (x,y) is of the form 
cx a {^ — r) d or cy a {-^ — r) d for positive integers a,m and some nonzero c, r, the first 
corresponding to the case where e lies below (d, d) and the second corresponding to where 
it lies above (d, d). Switching axes if necessary, assume that S e (x, y) = cx a (-^; — r) d . One 
can argue as in (7.3) — (7.8) to obtain a function of the form T(x,y) = S(x,y + a(x)) 
such that the bisectrix intersects N(T) at the vertex (d,d), but such that as in the last 
paragraph of the proof of Lemma 7.3 if (d, d) is the upper vertex of an edge e' of N(T) then 
T e /(1,7/) does not have a zero of order d (or greater). We also must have that T e >{— 1,7/) 
has no zero of order d or greater; for if it did by Lemma 7.1 we could write T e >(x, y) in the 
form cx a (^ — r) d , which would imply T e /(1, y) also has such a zero, a contradiction. 

There still remains the possibility that (d, d) is the lower vertex of a compact 
edge / of N(T) such that 7/(1,7/) has a zero of order order d or greater. By Lemma 7.1, 
if this happens Tf(x, y) is of the form cy a ( J j — r) d where g is an integer. Again using the 
argument from (7.3) onwards, this time reversing the roles of the x and y variables, there 
is a smooth coordinate change of the form (3 : (x, y) — ■> (x — y 9 h(x, y), y) such that if one 
denotes T o (3 by S' , then if (d, d) is the lower vertex of an edge /' of N(S') then Sf,(x, 1) 
and Sf,(x, —1) do not have any zeroes of order d or greater other than x = 0. This means 
S'f,(l,y) and S'*, (— 1,7/) also have no such zero. For if one of the functions did, S'j,(x,y) 
would have zeroes of order d along some curve y = sx n , which would imply either Sf,(x, 1) 
or Sf,(x, —1) had zeroes of order d or greater away from the origin, a contradiction. 

Furthermore, the slope of e' is of the form — for m' a positive integer and 
g > -^j. As a result, e! is an edge of N(S') containing (d, d) and the coordinate change 
(3 did not change any of the terms of T e /(x,y). Thus S' e ,(x,y) = T e >(x,y) and all zeroes 
of S' e ,(l,y) and S' e ,(—l,y) other than y = have order less than d. Hence we are in 
superadapted coordinates and the proof of Theorem 7.4 is complete. 
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